Interlayer Coupling in Ferromagnetic Semiconductor Superlattices 
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We develop a mean-field theory of carrier-induced ferromagnetism in diluted magnetic semiconduc- 
tors. Our approach represents an improvement over standard RKKY model allowing spatial inho- 
mogeneity of the system, free-carrier spin polarization, finite temperature, and free-carrier exchange 
and correlation to be accounted for self-consistently. As an example, we calculate the electronic 
structure of a Mn^Gai-^As/GaAs superlattice with alternating ferromagnetic and paramagnetic 
layers and demonstrate the possibility of semiconductor magnetoresistance systems with designed 
properties. 
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Semiconductors and ferromagnetic materials play a 
largely complementary role in current information pro- 
cessing and storage technologies. The possibility of ex- 
ploiting synergies between their properties has led to 
interestEl in fabricating hybrid systems. Among these, 
the most widely studied are ones in which ferromagnetic 
metals are patterned on the surface of a semiconductor E 
An attractive alternative has been presented by recent 
advancesflu in the fabrication and control of diluted mag- 
netic semiconductors which exhibit free-carrier induced 
ferromagnetism. In semiconductors, low intrinsic car- 
rier densities allow electrical transport properties to be 
grossly altered by changing doping profiles or gate volt- 
ages. In ferromagnets, it is long range order which al- 
lows weak magnetic fields to reorient large moments and 
produce significant changes in transport and other prop- 
erties. Systems with magnetically ordered free carriers 
in semiconductors open up a plethora of intriguing pos- 
sibilities for engineered electronic properties. Since these 
will generally involve spatial patterning of dopant and 
magnetic ion densities, their description requires a theory 
of free-carrier induced ferromagnetism in semiconductors 
appropriate to inhomogeneous systems. In this Letter we 
present such a theory and apply it to the case of layered 
structures w-hich mimic metallic giant-magnetoresistance 
multilayers. We conclude that strong magnetoresistance 
effects associated with half-metallic (fully polarized free- 
carrier band) states can occur in such systems. 

The systems in which we are interested are described 
by Hamiltonians of the form 



H = TL m + Hf - Jpd ^2 Si ■ srfin - Ri 



(1) 



i,I 



Here TC m is the Hamiltonian of localized spins Si, lo- 
cated at positions Ri and interacting with an external 
magnetic field and, in general, also with each other. TLf 
is the Hamiltonian of a semiconductopiree-carrier system 
described using an envelope functior£3 language and in- 
cluding interactions of carriers with a random disorder 



potential and with other carriers. We assume here that 
the semiconductor of interest has a single parabolic band; 
the formalism is readily generalized to multi-band situa- 
tions which will often be of practical interest. The last 
term on the right-hand side of Eq. (]l|) represents the ex- 
change interaction! 11 ] between Si and free-carrier spins 
Si which is responsible for the_ novel physics. Our the- 
ory is simplified by the neglectE] of dynamic correlations 
between a localized spin and both other localized spins 
and the free-carriers. We use a grand-canonical ensem- 
ble where the full system is in contact with a heat bath 
and the free carriers are in contact with a particle reser- 
voir. Since wCj-will ultimately use the local-spin-density 
approximationll3 (LSDA) to account for correlations in 
the inhomogeneous free-carrier system, it is convenient 
to formulate our mean-field-theory in terms of the Gibbs 
thermodynamic variational principle which asserts that 



0, < -k B T\n 



J2eM-{a\K\a)/k B T) 



(2) 



Here \a) ranges over all possible direct products of local- 
ized spin moment mi configurations and many-fermion 
free carrier states obtained by diagonalizing H in sub- 
spaces with fixed localized spin configuration. 



(a\H\a) = - ^{g^ B B)mi + E K [mi\ 



(3) 



where E K [mi] is a many-particle eigenvalue for 
free-carriers whose Hamiltonian includes the spin- 
dependent external potential term, — s z ^[g* fi B B + 
Jpd^2i Tni5(fi — Ri)] with g and g* denoting localized- 
spin and free-carrier g-factors respectively. 

The sum over the many-fermion free-carrier states 
can be performedJbrmally by appealing to spin-density- 
functional theoryEa which states that the free-carrier 
grand potential can be expressed in terms of the ground 
state number and spin densities. Assuming that the total 
spin quantum number of each localized moment is S, it 
follows that 
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n = ^nM(b I ) + n f [n h n l ]- 
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dfm(f) (4) 



where Qt [nj, nj is the internal free-carrier contribution, 
hi = [gn B B + J pd m(Ri)/2}/k B T, and 



ft (0) (6/) = -k B T\n [ exp(bii)] 



(5) 



is the thermodynamic potential of an isolated localized 
spin. In Eq. (|J), m(r) = n^(r) — n^(f) and the mean-field 
densities, n-\(r) and n^if), of the itinerant spins are to 
be adjusted so that the functional derivatives of Q with 
respect to both n-\(r) and ni(r) are identically equal to 
the chemical potential. 

In spin-density functional theory practical self- 
consistent field calculations for inhomogeneous inter- 
acting fcrmions are made possible by the Kohn-Sham 
separation of Qf into single-particle, electrostatic and 
exchange-correlation pieces. The development here is 
standard apart from the introduction of an additional 
spin-dependent effective potential because of the depen- 
dence of bi on n^(Ri) and ni(Ri). We find that the 
equilibrium free-carrier spin densities and the mean value 
of the localized spin quantum numbers on each site can 
be determined by solving the following equations self- 
consistently: i) a free-carrier single-particle Schrodinger 
equation with a spin-dependent potential, 



2m a 



(6) 



where 



V„(r) = v es (r) + v xc<a (r) - - (g*f-i B B + h pd (r)) ; 

ii) the Poisson equation for the electrostatic potential 
- 2 ' n(f) 



v es (r) — — I dr 1 



v ex t{r) 



(7) 



where 



n(f) = n T (r) +n l (f) , n a (f) = ^ f(ek,<r)\ipk,*(r)\ 2 , 

k 

f{^k,a) is the Fermi distribution function, -and v ex t{r) 
is the envelope function external potentialli!j including 
most importantly band edge and ionized impurity con- 
tributions; iii) the LSDA equation for the spin-dependent 
exchange-correlation potential 



d[ne 3:c (n T ,n x)] 
dn n 



\n a =n<r(f) ' 



(8) 



where e xc (ni,ni) is the exchange and correlation energy: 
per particle of a spatially uniform free carrier system;E£l 
iv) the mean-field equation for the exchange-coupling ef- 
fective Zeeman field 



h p d(r) = J P d^S(r- Ri)(mi) 



(9) 



where the mean-field localized-spin moment|— (mi) — 
SBs(biS) and B s (x) is the Brillouin function.lli 

For homogeneous systems with randomly distributed 
localized spins these equations can be solved analytically 
and represent an improvement over the mean-field theory 
RKKY interaction description of free-carrier induced fer- 
romagnetism. Our approach allows finite temperature, 
free-carrier exchange and correlation, and free-carrier 
spin-polarization effects to be conveniently accounted for. 
The Curie- Weiss temperature, obtained from Eqs. (^)- 
(§), is given by 



knT r . = 



cS(S+l) 



J 2 



-Xf(n,T) , 



(10) 



where c is the magnetic impurity density, and x/( n i T) is 
the temperature and density dependent free-carrier mag- 
netic susceptibility. In Fig. [I] we have plotted the fer- 
romagnetic transition temperature predicted by this ex- 
pression for p-type Mn^Gai^As with S = 5/2, J p d — 
0.15 eV nm 3 , c = 10 21 cm~ 3 and hole mass m* = 0.5m e , 
as a function of free-carrier Fermi wavevector fcp-. 
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FIG. 1. Curie- Weiss temperature, as a function of 
free-carrier Fermi wavevector, calculated including (solid line) 
and neglecting (dashed line) the exchange-correlation po- 
tential. The localized spin moment (upper inset) and the 
free-carrier relative spin polarization (lower inset) are plotted 
as a function of temperature for kF = 1.4 nm -1 . 

When exchange and correlation effects are neglected 
and x/(n, T) is replaced by its zero-temperature value, 
T c is proportionallo kp and agrees with the RKKY the- 
ory expression. QoEj Free-carrier exchange and correla- 
tion enhances T c by ~ 30 K in the range of free-carrier 
densities studied. The theory appears to be reasonably 
accurate when applied to the experimental systems we 
are interested in, those which have the largest critical 
temperaturesu (T c sa 110 K), providing confidence in its 
application to inhomogeneous systems. It is not adequate 
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at small Mn fraction (x < 0.02) where ferromagnetism is 
not observed, presumably because the free-carrier density 
is below its Mott limit. Experimentally, T c also decreases 
for x > 0.07, possibly because of soin-fluctuation effects 
neglected in the mean-field theoryO. The insets in Fig. |l| 
show that in the density-range of interest, the magneti- 
zation of the localized spins saturates more slowly than 
that of the free carriers as the temperature falls below 
the critical temperature. 
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FIG. 2. Interlayer exchange coupling (solid lines) as a func- 
tion of the Mn-undoped GaAs spacer thickness (a) and as 
a function of the average density of free-carriers (b). The 
dashed lines is the RKKY interaction coupling energy cal- 
culated with all magnetic impurities confined to planes sep- 
arated by the superlattice period. Results are plotted as a 
function of dimcnsionless parameter 2dkF- 

Once the band edge and ionized and magnetic im- 
purity profiles have been specified, Eqs. (||)-(^|) can 
be used to solve for the sytem's equilibrium proper- 
ties. As an illustration of our approach, we consider a 
Mn^Gax^As/GaAs superlattice with magnetic impuri- 
ties in alternate layers. We look for two different self- 
consistent solutions of Eqs. a ferromagnetic (F) 
one with parallel ordered moments in all Mn-doped re- 
gions, and a solution with an antiferromagnetic (AF) 
alignment of adjacent magnetic layers. The interlayer ex- 
change coupling E c , is defined as the difference in energy 
between AF and-F -states per area per Mn^Gai-^ As layer 
and is expecteefcj to oscillate with GaAs spacer width. In 
Fig. H we present numerical results for E c as a function of 
a dimcnsionless parameter 2dkF where kp is the Fermi 
wavevector corresponding to the average 3D density of 
free carriers in the superlattice with a period d. The 
system we consider has 2 nm thick Mn^Gai-^As layers 
with c = 10 21 cm -3 and homogeneously distributed ion- 
ized impurities that neutralize free-carrier charge. Oscil- 
lations in the self-consistent (solid lines) E c , are quality, 
tively consistent with simple RKKY model estimates ;E£l 
differences exist primarily because of the proximity in- 
duced spin-polarization in the nominally paramagnetic 



GaAs layers mentioned below. The amplitude of oscillsu. 
tions in E c is ~ 10 x smaller than in metallic systemscll 
measured in absolute units and ~ 10 x larger if energy 
is measured relative to the Fermi energy of free carriers. 
In order to achieve substantial exchange coupling in ex- 
perimental systems it will be important to limit disorder 
scattering in the GaAs layers. 
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FIG. 3. Self-consistent results plotted within the unit cell 
of the Mn^Gai-zAs/GaAs superlattice. From top to bottom: 
i) effective potentials from Eq. ^ for up (solid lines) and down 
(dashed lines) spin free-carriers (the dotted line here indicates 
the chemical potential); ii) mean-field localized spin moment; 
iii) free-carrier relative spin-polarization. 

Our mean-field calculation also yields information on 
the localized and free-carrier spin magnetization densities 
and on the spin-split bands of the free-carrier system. 
This information provides a starting point for building 
a theory of electronic transport. Numerical results for 
the above Mn^Gai-^As superlattice with d — 3.5 nm 
and kp = 1.4 nm -1 (corresponding to average density 
10 20 cm -3 ) are summarized in Figs. || and |j. In the 
top panels of Fig. || the chemical potential and the ef- 
fective potentials for spin-up and spin-down free carriers 
(see Eq. (^)) are plotted as a function of z over a AF 
configuration unit cell. All energies here and below are 
measured from the spatial average of the electrostatic 
potential v es (z). The potentials V a have similar shapes 
in F and AF cases, except for the reversed order for up 
and down spins in the right Mn-doped layer. Note that in 
this example, confinement of carriers in the magnetic lay- 
ers is due entirely to the exchange potential produced by 
magnetic order. The localized (middle panels) and itiner- 
ant (lower panels) spin systems reach 100 % polarization 
in the Mn-doped layers at the temperature (T = 0.1T C ) 
and carrier density for which these calculations were per- 
formed. The itinerant system spin-polarization is large 
in the layers free of magnetic impurities, especially so in 
the F-state case. 

Fig. |] shows occupied minibands in the superlattice 
Brillouin zone. In the F-state spin-up and spin-down 
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minibands are split by about 0.25 eV. There is no spin- 
splitting in the AF-state since the effective potentials V\ 
and V[ differ by a rigid shift in the z direction. The mini- 
band dispersion is much weaker in the AF case because 
the barriers separating two adjacent minima in the effec- 
tive potential are twice as thick and high as in the F case. 
Since the conductance is approximately proportional to 
the square of the largest miniband width in either coher- 
ent or incoherent transport limits, the minibands can be 
used to estimate the size of the current-perpendicular- 
to-plane (CPP) magnetoresistance effect. For the case 
illustrated, the AF state CPP conductance will be three 
orders of magnitude smaller than the F state CPP con- 
ductance. The large difference is expected since the bulk 
Mn^Gai-ajAs bands are half- metallic for these param- 
eters. In general we predict strong CPP magnetoresis- 
tance in Mn^Gai-^As/GaAs multilayer systems if the 
AF state can be realized. 
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FIG. 4. Partially occupied energy minibands in F-state 
for spin-up (solid line) and spin-down (dashed fine). In the 
AF-state both spins (dotted line) have the same minibands. 
The chemical potential is 0.053 eV and a— 7 nm is the unit 
ceif fength. 

The GaAs/Mn x Gai_a;As superlattice is but one ex- 
ample of a semiconductor nanostructure to which our 
mean field theory can be applied. Tremendous flexibil- 
ity is possible even within systems containing only GaAs 
based materials. These can act like a normal (N) con- 
ductor, like an itinerant ferromagnet (F) when a fraction 
of the Ga atoms is replaced by Mn, and like an insulator 
(I) when a large enough fraction of Ga in the crystal is 
replaced by Al. In GaAs-based N/F/I heterostructures 
all components are lattice matched and have relatively 
simple band structures. Geometries of interest include 
F-N junctions, N-F-N spin-filters, F-N-F spin- valves, and 
F-I-F magnetic tunnel junctions. Since the materials are 
semiconductors rather than metals, external bias voltages 
can have a strong influence on both charge and magneti- 
zation density profiles. Self-consistency is hence essential 
in modeling these systems. 
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